Comparison of Nuclear Suppression Effects on Meson Production at High px and pi 
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The medium effect on the pion distribution at high pr in AA collisions is compared to that of the 
pion distribution at high pl in pA collisions. Both the suppression of the spectra and the energy 
losses of the measured pions are studied. Although the medium effect on pr is larger than on pl, 
the difference is found surprisingly to be not as big as one would naively expect. 
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I. INTRODUCTION 

When partons traverse nuclear medium, whether dense 
or not, they lose momenta by scattering and gluon radia- 
tion. At high transverse momentum px the effect can be 
calculated in perturbative QCD, although its reliability 
is not expected for pr < 10 GeV/c. At low pr, but high 
longitudinal momentum p^, the effect cannot be calcu- 
lated in pQCD, although it is known that partons suffer 
momenta losses also in the beam direction. Experimen- 
tally, it is the momenta of the produced hadrons that 
are measured. How they are related to the underlying 
parton px and pl distributions is still controversial. But 
even at the phenomenological level it is unknown what 
the relationship is between the properties of momentum 
degradation in the transverse and longitudinal directions. 
We attempt to shed some light on that relationship in this 
paper. 

There are currently good data on high-j»T ^° produced 
at the Relativistic Heavy-Ion Collider (RHIC) up to 
pr = 8 GcV/c at various centralities dj. One can there- 
fore deduce the suppression factor from the pr distribu- 
tions as a function of px and centrality. There are no 
comparable data on the pl distributions from RHIC. At 
lower energies the data of NA49 at the SPS provide a 
good description of the effect of baryon stopping in pA 
collision What we need for comparison with the pr 
of 7T° is the pl distribution of the produced pions, for 
which no data are yet available. However, we do have a 
model calculation of the pion distribution that contains 
the degradation effect extracted from the observed pro- 
ton distribution that will be our input in our study 
of the suppression factor in the pi, distribution. 

The two features that we shall compare are very differ- 
ent: the degradation of pr in Au + Au collisions at RHIC 
and the degradation of pl in p + Pb collisions at SPS. 
In the absence of any information in the literature on 
the quantitative or qualitative difference between the two 
types of degradation effects of the nuclear medium, even 
a crude estimate of the suppression properties described 
in a common language would be illuminating. Without a 
study of the type proposed here, one does not even know 
whether the strengths of suppression are within the same 
order of magnitude, especially since the medium in AA 



collision at RHIC is dense and hot, while the medium in 
pA collision at SPS is uncompressed and cold. 

There is a theoretical issue that is of interest to discuss 
here. For some time there has been a school of thought 
that all hadrons are produced by the fragmentations of 
partons, not only in the transverse direction in the form 
of jets (which is generally accepted), but also in the lon- 
gitudinal direction in the form of breaking of strings (as 
in the dual parton model The effect of the nuclear 
medium has been represented by the modification of the 
fragmentation function |5j in transverse direction, but 
that of the longitudinal direction is not known. However, 
fragmentation of partons is not the only way to produce 
hadrons. Recent investigations have shown that quark 
recombination can be important in the high px problem 
0-0- i n addition to its relevance originally proposed for 
the high p^ problem [slllOj. Since the multiparton distri- 
butions needed for recombination are drastically different 
in the transverse and longitudinal directions, the effect of 
the nuclear medium is considerably more complicated. In 
this paper we can avoid dealing directly with those com- 
plications by putting the emphasis on the phenomenology 
of the hadrons produced. For the yj problem we shall 
use the scaling form of the data [Tj, while for the pl 
problem the calculated x p distributions that follow from 
the pA data will be used. Our task is made easier by not 
deriving the pt and pl distributions, but by focusing on 
the centrality dependences of those distributions. 

Suppression of the meson distribution is like jet 
quenching at the parton level. In addition to the study 
of suppression, we shall also consider energy loss, which 
is another way of quantifying the medium effect. It cor- 
responds to a shift in px or p^ that is necessary for the 
inclusive cross section in medium to be equivalent to a 
reference cross section with minimal medium effect. Wc 
shall find interesting results in the shift that are very dif- 
ferent from the prediction of pQCD. That difference is 
much larger than the difference between the effects on 
the transverse and longitudinal motions of the produced 
mesons. 

Since our comparison is between pr in AA collisions 
and in pA collisions, they are two steps removed from 
each other. When good pr, data on identified pions in 
AA collisions become available, they shall then serve as 
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the intermediate station to make possible two one-step 
comparisons: (1) between pt and pi, in A A collisions, 
and (2) between pi, in AA andp^ vnpA collisions. What 
we do here therefore sets the stage for that work to come. 



II. SUPPRESSION OF THE PION p T 
DISTRIBUTION IN AA COLLISIONS 

For the pr distribution of pions in nuclear collisions we 
use the PHENIX data on tt° production at midrapidity 
with px extending to as high as 8 GeV/c l]. A conve- 
nient scaling form for that distribution has been found 
that summarizes the dependence on energy and central- 
ity in terms of a simple analytical formula 0, ^| . The 
quantity that we want to study is the suppression factor 
S, which is the ratio of the normalized px distribution 
P: 



S(x T ,N)=P(x T ,N)/P(x T ,2) 
where xt is the scaled px variable 
x T = p T /K , 



(1) 



(2) 



and N is the abbreviated notation for the number of par- 
ticipants iVp art . The scale Kq is set at 10 GeV/c for con- 
venience; it is trivial to move it higher when higher px 
data become available. If x^ 1 dN„ / dxT denotes the xt 
distribution of produced 7r°, averaged over midrapidity 
and over all azimuthal angle cj>, then P(xt, N) is defined 
by 



P(x T ,N) 



1 dN„ 

xt dxr 



dxi 



dNi, 
dxr 



(3) 



Instead of determining P(xt, N) directly from the data 
for every N at any given s, it is simpler to derive it from 
the scaling function $>(z) that is an excellent fit of all 
high-energy data for all centralities 0, ^]. Since the 
details of the relationship between P(xt,N) and $(z) 
are given in Refs. 0, 0, there is no need for us to 
repeat them here. For the convenience of the reader, 
they are summarized in the Appendix. 

From the analytical expression for S(xt,N) obtained 
by use of Eq. (|60[) in Eq. we can examine the N 
dependence for fixed values of xt by plotting \iiS(xt, N) 
vs N for some sample values of xt, as shown in Fig. 1. 
Approximating the nearly linear behaviors in Fig. 1 by 
straight lines, we obtain 



ItiS(xt, N) = ao(xr) — cli{xt)N , 



(4) 



where oo(xt) and ai(xr) are shown in Fig. 2. The solid 
lines in that figure are fits, using the parametrization 

a (x T ) = 0.031 [1 - exp(0.5 - 12x T )] , (5) 



ai (x T ) = 0.0043 [1 - exp(0.2 - 5x T )] . (6) 
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FIG. 1: The dependences of the suppression factor S(xt,N) 
on the number of participants iV for various values of xt- 




FIG. 2: The coefficients clo(xt) and ai(xr), where o,o(xt) is 
plotted with reference to the scale on the left, while o,i(xt) 
refers to the scale on the right. 



A general statement that can be made about IiiS(xt, N) 
is that it is approximately linear in N, and that the pa- 
rameters of the linear fits are roughly independent of xt 
when xt _ 0.4. 

Since the number of participants N is associated only 
with AA collisions, we must express it in terms of some 
measure of nuclear length in order to be able to compare 
the suppression factors in A A and pA collisions. The 
dependence of N on the impact parameter b is known 
|13| . At a fixed b, the lens-shaped overlap region in the 
transverse plane of two colliding nuclei of the same radius 
Ra has a minimum distance between the center of the 
overlap and the edge of either nuclei (assumed to have a 
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sharp boundary) 

imin = R A - 6/2 • (7) 

The maximum distance between the center of the overlap 
to the edge of both nuclei is 

L max = {R 2 A ~b 2 /A) 1/2 . (8) 

The distance that a parton would travel in the nuclear 
medium at midrapidity in the transverse plane can be as 
large as 2L max and as small as 0, depending on where 
the parton starts and in which direction. We shall set 
the average distance L traversed to be 

L(b) = L min = R A - 6/2 , (9) 

which is an approximate average over all azimuthal angles 
and origins of parton paths. Any more detailed geomet- 
rical averaging is pointless, since the nuclear density of 
the overlap region is not uniform, and the nonuniformity 
depends on b. Knowing N(b) and L(b) enables us to plot 
N vs L, as shown by the solid line in Fig. 3. We shall 
approximate N(L) by a straight line 

N(L) = -35.8 + 62. 2L (10) 

(with L in units of fm), shown by the dashed line in Fig. 
3. In view of the nonuniformity of the nuclear density 
traversed by a parton in the overlap regions, we believe 
that an approximation of N(L) by a linear dependence in 
Eq. QlOfl is good enough to represent the path length that 
enters into the description of the momentum degradation 
effect. 

400 1 1 1 1 > — i 




FIG. 3: The parametric dependence of N(b) on L(b), where 
the solid line follows from the use of Eq. for L(b), while 
the dashed line is a straightline fit of the solid line. 

Substituting Eq. (|TU|) into Eq. Q we obtain 

lnS(x T ,L) = A (x T ) - A(x T )L , (11) 



where in Ao(xt) the contribution from the o,q{xt) term 
in Eq. (J1J is negligible. Thus Eq. (|llf> may be rewritten 
as 

]nS(x T ,L) = -A(x T )(L - 0.71) , (12) 

where 

A(x T ) = 0.27 [1 - exp(0.2 - bx T )] , (13) 

which is nearly constant for xt > 0.5, as shown by 
the solid line in Fig. 4. This exponential dependence 
of S(xt, L) on L is in a familiar form for nuclear attenu- 
ation. We remark that although the length L is the ap- 
proximate average distance in the nuclear medium that a 
parton traverses, no parton dynamics has been assumed 
in the derivation of S(xt, L), which is extracted from the 
data on ir° production without dynamical modeling. The 
exponential form of S(xt,L) is similar to the Gerschel- 
Hiifher formula for the J/ip suppression , except that 
the latter refers to a quantity integrated over all px, and 
is related to the dissociation of J/ i/j in the medium. The 
quarks and antiquarks that form the produced ir° in our 
expression for S(xt,L) are the results of gluon radia- 
tion and gluon conversion processes, most of which are 
not calculable in pQCD, especially near the end of the 
evolution process. 
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FIG. 4: The solid line is for A(xt) and the dashed line is for 
X(xl)- The horizontal axis can be for either xt or xl- 



III. SUPPRESSION OF THE PION p L 
DISTRIBUTION IN P A COLLISIONS 

Since no distributions of identified particles in the 
fragmentation regions of AA collisions at RHIC are avail- 
able, we consider the processes at the SPS energies. At 
lower energies the projectile fragmentation region can 
contain particles arising from the fragmentation of the 
target, and vice- versa. It is therefore important to con- 
sider reactions that are free of such 'spill-over' particles. 
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The preliminary data of NA49 on pA collisions at SPS 
have been analyzed to provide the produced p — p distri- 
bution of the projectile fragmentation only Q- That is 
accomplished by determining the pl distribution 

F P {p-p) = (p-p) P - 2 [(P-P)n+ + (p-p)*-] , (14) 

where (p—p)h denotes the distribution of h+A — > (p— p) + 
X, since the quantity inside the square brackets contains 
no beam fragments by charge conjugation symmetry. A 
similar distribution for pions, F p (ir + — fl" - ), unfortunately 
does not exist, which is what we need for comparison with 
the result obtained in the preceding section. 

Although there are no experimental data for pion pro- 
duction in the proton fragmentation region free of target 
fragments, theoretical results on such distributions are 
available that are based on the centrality dependence of 
the experimental data of F p (p— p). In Ref. [3j the recom- 
bination model is used to relate the data to the nuclear 
degradation effect on the pl of the produced p — p, which 
in turn is then used to predict the pl distributions of the 
produced pions. Since the nuclear suppression factor on 
the produced pions is directly related to the experimen- 
tal suppression of the produced p — p, we shall use the 
result in Q for the pion distribution in p + Pb collisions. 

Let us use H(xl, v) to denote the pion inclusive distri- 
bution in the scaled longitudinal momentum, integrated 
over pt, i.e., 



H(x L ,v) = x L 



dxr. 



(15) yields 



where xl = %Pl/ \fs and v is the average number of col- 
lisions with nucleons in the target nucleus. From Fig. 9 
in Ref. Q one sees that H(xl, v) is essentially exponen- 
tial in ccl for various values of v. We use the following 
parameterization for p + Pb — > tt + + X 



H(x L , v) = exp [-fto(f) - h\(v)xL\ , (16) 



where 



V3.1) = 0.92, h a (6.3) = 1.61 , 

/ii(3.1) = 5.89, h x (6.3) = 6.57 . (17) 

These parameters provide a good fit for 0.2 < xl < 0.6. 
For xl > 0.6, the slopes are slightly higher, but the ratio 
of H(xl, v) at the two values of v is about the same. 

Since the data are for the two values of v only, our 
definition of the suppression factor is 



S (x L ;v 2 ,vi) = H (x l ,l> 2 ) /H(xl,vi) 



(18) 



with V\ — 3.1 and v 2 = 6.3. We now make the assump- 
tion that ho(D) and h\{v) are linear functions of v, partly 
because the average (xl) can be shown to decrease expo- 
nentially with v [3] , and partly because all nuclear damp- 
ing effects are empirically dependent on the path length 



in exponential form. With that assumption we can ex- 
press S (xl; ^,1) for any v relative to v = 1 as 

lnS (x L ;v, 1) = - Mv) + Q x {v)x L \ (19) 

where Oq{v) and 9±(D) are the linear functions of v 

6i{y) = hi[y)-hi(l) =6' i {v-l) , i = 0, 1 . (20) 

The slopes are 

= hj{D 2 ) ~ hjjpx) 

From the values given in Eq. I|17|l , we find that 0' and 6[ 
are very nearly equal. We denote them collectively by 

9' = 0.21 (22) 

and Eq. (fTU)) becomes 

\nS(x L ;DA) = -9'(l + x L )(D-l) . (23) 

The exponential dependence of S on v is now explicit. 
The specific xl dependence of the decay coefficient is a 
result of the distribution in Eq. H16[) calculated in jjj. 

If L denotes the path length in a pA collision, then the 
usual expression for v in terms of L is 

47T A . , 

u = a pp pL, P =Y~R\' ^ ' 

Using Ra — 1.2A 1 / 3 fm, and a pp = 35 mb, we have 

v = 0.48L , (25) 
where L is in units of fm. Substituting this in Eq. 1231) 



]nS(x L ,L) = -\(x L ) (£-2.1) , (26) 



where 



\(x L ) = 0.1 (1 + x L ) . (27) 

Equations (|26|l and l|27|) for the suppression of the pl 
distribution in pA collisions are the counterparts of Eqs. 
(|12[) and (|13fl for the suppression of the distribution 
in AA collisions. Instead of A(xt) that rises rapidly at 
small xt, but saturates to a nearly constant value for 
xt > 0.5, as shown in Fig. 4, we now have X(xl) which 
is linearly rising, shown by the dashed line in Fig. 4. Since 
Eq. H16(l is not reliable for xl < 0.2, we do not show that 
portion of \{xl) in Fig. 4. For numerical comparison we 
can consider xt — xl at two values, 0.6 and 0.8: 



A(0.6) = 0.25 . 
A(0.8) = 0.26 



A(0.6) = 0.16 . 
A(0.8) = 0.18 



(28) 



Evidently, A(xt) is larger than \{xl) when xt = xl > 
0.5, but not by much. Of course, there is no cogent reason 
to compare them at xt ~ xl since the scale Kq in the 
definition of xt is arbitrary, while xl is Feynman xf 
and has scaling property. However, with A(xt) being 
roughly constant for xt > 0.5, it does not matter what 
xt is exactly. The relative values of A(xt) and \(xl) 
shown in Eq. (|28|l give us some indication of how they 
differ. 
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IV. ENERGY LOSSES IN p T AND p L 

Another way to quantify the medium effect is in terms 
of energy loss. Since our aim in this paper is to stay at 
the level of observable quantities, we cannot descend to 
the parton level where the concept of energy loss makes 
sense as one can study the evolution of a parton in its 
trajectory through the nuclear or quark medium. A pro- 
duced meson does not itself traverse that medium, since 
it is only at the end of the evolution of the parton system 
that it is formed. Nevertheless, the notion of energy loss 
can be expressed in terms of a shift in xt in comparing 
the meson inclusive distributions at two different values 
of N. To be specific, let us set the reference value of N 
at N = 2, which is not exactly pp collision, but is low 
enough to represent minimal nuclear effect. Let us then 
define the shift X in xt by 



P(x T ,N) = P(x T +X,2) 



(29) 



where P(xt,N) is the normalized pion distribution de- 
fined in Eq. |JSJ. Thus X measures the degradation of 
the pion xt in changing iV par t from 2 to N. 

Using Eq. (|50|) where (xt) is calculable by use of Eq. 
(JSBJ, we can solve Eq. (|29|l and determine X in terms of 
xt and N. At fixed N, the dependence of X on xt is 
nearly perfectly linear. Paramctrisizing it as 



X(xt,N)=x (N)+Z(N) x T 



(30) 



we find, for N = 100,200, and 350, the result 
x (N) = -0.0047,-0.0023, and -0.0039, and £(N) = 
0.0432,0.0916, and 0.1741, respectively. Since x is neg- 
ligible except at very small Xt, we can regard £(N) as 
the fractional shift, or more precisely as 



£(JV) = dX(x T ,N)/dx T , 



(31) 



whose dependence on TV is shown in Fig. 5. Disregarding 
the point at N = 2, where by definition £(2) = 0, the 
three points at N = 100, 200, and 350 can be fitted by 
a straight line, as shown by the solid line in Fig. 5. Its 
parametrization is 



£(jV) = -0.011 + 5.26 x 10~ 4 N 



(32) 



for N ~ 50. Thus for AA collisions the fractional energy 
loss or fractional shift in pt is independent of pt and 
depends linearly on N with a coefficient 



(,' N = 5.26 x 10" 



(33) 



This result can be obtained quickly, but only approxi- 
mately, in two ways. First, if one ignores the dependence 
on N of the denominator in Eq. H53f) . an approximation 
that amounts to setting k(s,N) to be a constant (an er- 
ror < 20%) compared to the many orders of magnitude of 
variation of $(z), then the condition of Eq. (|29|) is equiv- 
alent to identifying its two sides with the same function 
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FIG. 5: The N dependence of the fractional shift £. The 
solid line is a linear fit of the upper three points. 



3>(z) with z evaluated at N pax t = N and 2, i.e., 
PT (2)-p T {N) K(2)-K(N) 



pr(N) K(N) 
6.36 x 10~ 4 (N — 2) 



K(N) 

Setting K(N) « K{2) = 1.226 results in 



N 



5.19 x 10" 



(34) 



(35) 



which is only slightly less than the more accurate value in 
Eq. (|33|l . The other approximate way of determining the 
fractional shift is to use Eq. 160fl and ignore the N de- 
pendence of (xt) in that equation so that the condition 
in Eq. I|29|l becomes an identification of the u variable in 
&(u) at two values of N, i.e., 



u{N) 



u(2) 



x T + X 



Since (xt) n depends on TV in a known way [Tl| 

(x T ) N = (x T ) No exp[-0(N-N o )] , 



(36) 



(37) 



where = 5.542 x 10~ 4 (denoted by A in [TJ), its use m 
Eq. I|36|) results in 



£' N « (3 = 5.54 x 10~ 4 . 



(38) 



Note that the smallness of (3 in Eq. H37JI roughly justifies 
the treatment of (xt) as a constant in Eq. I|60|l in the first 
place. The value of ^ in Eq. I|38(l is only slightly larger 
than that in Eq. I)33|l . Thus the two approximate meth- 
ods yield results that bracket the correct value closely, 
and illustrate the crucial role that the scaling variables z 
and u play. 

It is of interest to note that in pQCD the shift in 
Pt of the vacuum spectrum necessary to effect the in- 
medium spectrum is proportional to yJp T |l5j | . It means 
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that the fractional shift decreases with px, whereas our 
phenomenological result indicates that it is independent 
of pt- However, since pQCD is reliable only for px > 10 
Gev/c, while our analysis is based on data at pr < 10 
Gev/c, there is as yet no direct conflict. Nevertheless, 
the disagreement in the pt dependences is worth bearing 
in mind. 

We now consider the energy loss of pl of pions pro- 
duced in pA collisions. We define the shift X by referring 
the inclusive cross section at v to that at v = 1, i.e., 

H{x L ,v) = H(x L +X,l) . (39) 

The use of Eq. (HHJ) and (HOJ leads to 

X(x L ,D)={(p)(x L + l), (40) 

where 

e(*) = &(*-i), (41) 

C P = 6'/hx{l) = 0.0385 . (42) 

To compare with the result from AA collisions, let us 
convert both N and v to the average path length L. Us- 
ing Eq. Jinj in l|3^]l. we get 

£,t(L) — —0.03 + 0.0327 L . (43) 

Using Eq. in |gT), we get 

£ L {L) = -0.0385 + 0.0185 L . (44) 

In both cases the fractional shifts depend linearly on L 
with the coefficients 

£t = 0.0327 , g L = 0.0185 . (45) 

The ratio of these two coefficients is very nearly the same 
as the ratio of A(xt) to A(xl) in the region xt,l ~ 0.5 
[cf. Fig. 4]. Thus the study of energy loss and that 
of suppression give comparable results on the effects of 
the nuclear medium. The advantage of using Eq. I|45|l 
for comparing the nuclear effects on the transverse and 
longitudinal motions is that the fractional shift is inde- 
pendent of the scale Kq used in the definition of xt in 
Eq. (J2J. Thus the numerical values of £'n>£'v>£t an d £l 
are simple quantitative results of this investigation that 
can be reexamined in the future when experimental data 
at different energies for different colliding nuclei become 
available. 



V. CONCLUSION 

We have studied the suppression effect of the medium 
and the energy losses of the produced particles. The 
former is a comparison of the inclusive distributions in- 
medium versus minimal-medium at the same xj- or x^. 
The latter is the shift in xt or i/, necessary for the two 



distributions at different N to be equivalent. Both stud- 
ies yield qualitatively the same level of effect. In the 
following we shall use the suppression effect to repre- 
sent both in our discussion of the differences between 
the transverse and longitudinal effects. 

The primary remark to make is that the values of 
A(xt) and X(xl) given in Eq. i|28|) are amazingly close. 
There are many arguments one can give to suggest that 
A(xt) and \(xl) should not be similar in magnitude, and 
few are available to explain that they are even within the 
same order of magnitude. Let us present some of them. 

The main difference between the two suppression ef- 
fects is that one refers to transverse, the other longitudi- 
nal motion. For transverse momenta of partons caused 
by hard collisions, at least there is pQCD to describe 
some aspects of the dynamics, although not reliably for 
Pt < 8 GeV/c. For the longitudinal momenta of the par- 
ticles detected, there is no basic theory to describe their 
behavior in terms of partons without some substantial 
use of models. Some properties of the pr degradation can 
be calculated, but a recent result on the nuclear modi- 
fication factor does not reproduce even the trend of the 
Pt dependence [lfj. Nevertheless, much more work has 
been done in applying pQCD to the high px problem 
than to the high pi, problem. Since no large momentum 
transfers are involved at high p l , one would naively not 
expect the nuclear suppression effect to be of the same 
nature as at high px- Yet we find A(xt) and X(xl) to be 
comparable. 

The media of the two problems are also different. In 
AA collisions at RHIC one expects the nuclear medium 
to be dense and hot, if not a quark-gluon plasma. In 
pA collisions at SPS the medium that the projectile tra- 
verses is the normal uncompressed nucleus. One would 
therefore expect the effects of the media on momentum 
degradation to be very different, yet they are not. 

The estimates of the average path length L involve dif- 
ferent approximations for the two cases, and it is difficult 
to assess the effects of those approximations. The only 
way that px in A A collisions can be compared to pl in pA 
collisions is in the common language of exponential decay 
in terms of a path length L. To improve on this problem, 
we have to gain more information from experiments at 
intermediate steps by measuring p^ in AB collisions with 
various nuclear sizes of A. 

We can think of one reason that could possibly explain 
the closeness of A(xt) to X(xl). The inclusive cross sec- 
tions that we examine for the calculation of the suppres- 
sion factor S are for pions, not nucleons or other baryons. 
Whereas leading baryons are strongly related to the va- 
lence quarks in the projectile, the pions are more asso- 
ciated with the gluons, which undergo conversion to qq 
pair before hadronization. The depletion of gluons in 
the nuclear medium can lead to pion suppression in any 
direction. Evidence for gluon depletion even in pA colli- 
sions can be found in the suppression of J/tp production 
at large xl fl7Lll£j. Our result on the closeness of A(xt) 
and \(xl) may well suggest that gluon depletion is the 
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main mechanism for the suppression of both px and pl 
distributions. 

In quantitative terms it must be recognized that A(xt) 
is undisputatively larger than \{xl) and therefore pro- 
vides some comfort that the suppression effect is en- 
hanced when the medium is denser and hotter. What is 
unexpected is that it is not an order of magnitude larger. 
In order to fully understand the suppression problem we 
need a whole set of experiments that measure identified 
pions at high p? and pi, for all combinations of nuclear 
sizes in AB collisions at high energy. It is important to 
discover what is universal in the nuclear effects on the 
produced particles, and what is not. The finding in this 
paper constitutes an interesting and intriguing beginning 
in that direction. 



Acknowledgment 

We acknowledge helpful discussions with H. Huang and 
A. Tai at the early stage of this work. This work was 
supported, in part, by the U. S. Department of Energy 
under Grant No. DE-FG03-96ER40972. 



Appendix 

In this Appendix we summarize the basic formulas that 
relate P(xt,N) to two scaling functions found in [Tl| . 
The first scaling function is 



$(z) = 1200 (V + 2)- 4 8 (1 + 25 e 



-4.5z\ 



(46) 



which describes the px distributions at all N (number of 
participants) and all y/s in terms of one scaling variable 

z=p T /K(s,N), (47) 

where K(s,N) = K(s)K(N), with 0G3 

K{s) = 0.69 + 1.55 x 10~ 3 yi , (48) 



K(N) = 1.226- 6.36 x 10~ 4 A , 



(49) 



\/~s being in units of GeV and K(s, N) in units of GeV/c. 
$(z) is related to the xt distribution by 



= A(N)k 2 (s, N) 



1 dN„ 



(50) 



XT dxT 

where rapidity density is implied, and 

A(N) = 530A C (A)" - 9 , N C (N) = 0.44A 1 ' 33 , (51) 
and 

k(s,N) =K(s,N)/K , (52) 

Kq being an arbitrary scale, fixed at 10 GeV/c for the 
definition of Xt in Eq. J5J . It is a phenomenological fact 



that the combination of separate factors on the RHS of 
Eq. H50|) that individually depend on ij, N and s results 
in a universal function that depends explicitly on the one 
variable z only. 

Using Eq. JSDJ in Eq. © yields 



P(xt,N) = $(xt,N) I J dx T x T <f>(xT,N) (53) 
where $(z) is expressed in terms of Xt and N through 



z = x T /k(s,N) 



(54) 



Dependences of P and $ on s will not be shown explicitly. 
With Eq. (jSSJ, the distribution P{xt,N) can thus be 
analytically calculated. However, instead of performing 
the integration in the denominator for every N, there is 
an even simpler relationship that makes use of another 
scaling function. 

It is found in Ref. that there exists another scaling 
function 



*(u) = $(z(u)) / duu$(z(u)) , 



where 



Pt 



(p'l 



XT z 

W) = ~) 



(55) 



(56) 



The new scaling variable u endows ty(u) with a prop- 
erty that is analogous to the Koba-Nielsen-Olesen (KNO) 
scaling [ljj. The average (z), defined by 



J dzz 2 $(z) 



(57) 



J dzz§(z) ' 

is a constant (z) = 0.414, and is related to (xt) by 

(z) = {x T )/k(s,N) (58) 
due to Eq. where (xt) is defined by 



. dN„ I f dN* 

[Xt) — / dXTXT—, — / / dxT- 



(59) 



dxT / J dxT 

A comparison between Eqs. I|53|) and (|55|l yields 

= (xt) 2 P(xt,N) . (60) 

It is the combination of Eqs. (j56(l and (f f>Cj|) that has led us 
to regard ^(u) as a KNO-type scaling. The advantage of 
dealing with u instead of z is that k(s, N) is not explicitly 
involved in relating u to the observable pr- Evaluating 
Eq. we have H3 



= 2.1 x 10 4 (u 2 + 11.65)" 4 - 8 (1 + 25e 



-1.864m\ 



(61) 



Using Eqs. I|58|l and (|60|l . we now have an algebraic for- 
mula for P(xt, N), which can be used directly in Eq. (JJJ 
for the suppression factor S(xt-, N). 
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